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THEEE-DIMERSIOHAL SUPERSORIC ROZZLES ARD IREETS OF 
ARBITRARY EXIT CROSS SECTIOR 
By Jolm C. EwHjrd and Stephen H. Maslen 


SUMMARY 

A method is presented for obtaining three-dimensional imsymmetric 
supersonic nozzles and inlets from knoim axl symmetric flows. Stream- 
lines hounding the desired exit shape are traced throu^ the known basic 
flow solution to give the required unsymmetric wall contours. Several 
examples are given. 


IRTRODUCnOR 

Because of the great complexity of the method of characteristics 
in three dimensions (see, for example, references 1 and 2) , the accurate 
theoretical design of general three-dimensional supersonic nozzles and 
diffusers has not been feasible. The designer is thus limited to two- 
dimensional eind axl symmetric flows. In many cases this may he a serious 
disadvantage. In the case, particularly, of a hypersonic wind tunnel, 
it would he very helpful if a three-dimensional expansion coiild he 
used to avoid the usual thin slit which would appear as the throat of 
such a nozzle if it were two dimensional. Etirthermore, such a nozzle 
ml^t have better secondary flow characteristics than the corresponding 
t^TO-dimensional nozzle in which transverse pressure gradients may cause 
cross flow in, and local huild-up of, the boundary layer. Bowever, the 
usual three-dimensional nozzle is axially symmetric and thus precludes 
the use of flat schlleren windows. In the design of supersonic inlets, 
considerations of space, as well as of permissible flow turning, makp 
the use of shapes other than two-dimensional or axisymmetric ones 
desirable. 

A simple method for obtaining three-dimensional nozzles from 
axi symmetrical ones has been developed at the HAOA Lewis labooratoiy 
and is presented herein. The procedure is exact within the limitations 
of the method of characteristics. 
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METHOD 

In any inviscid fluid :fLow, the streamsheets, hy definition, form, 
surfaces across -which there is no flow and hence may he replaced hy 
solid boundaries. This fact -will herein he applied -to find unsymmetrlc 
nozzles emd inlets from axially symmetric flows. Two salient features 
of an axL symmetric flow are that the streamlines lie in planes throu^ 
the streamwise axis and that the flow in any one such plane is the 
same as that in any other. 

In order to find the nozzle of a desired shape, the axL symmetric 
nozzle ha-ving the desired length and Ifech number characteristics is 
computed and then, choosing the desired cross section at some station, 
for example, the exit, the stioamlines which pass throu^ the periphery 
of that section are fo-und. streanOlnes can he determined either 

hy constructing them from the kno-wn local flow directions or hy the 
better method of integrating the axlsymme-fcric stream function (appendix). 
The streamsheet formed hy these stieamlines will constitute -fche -walls 
of the desired inlet or nozzle. 

If the initial and final flo-ws are uniform and axial, it can 
readily he sho-wn -that the cross-section shapes at -the inlet and exit 
•will he similar. For any arbitrary streamline, the ratio of the radii 
at the initieil and final stations is given hy the one-dimensional 
continuity eq.mtion for Isentropic flow: 

r+1 



where M and r are Mach number and radius, respectively, and the sub- 
scripts i and f refer -to the initial and final s-tations. This 
ratio depends only on the initial and final Mach numbers and hence is 
Independent of the initial streamline. Likewise each streamline con- 
tln-ues in its own axial plane thro-ughout the nozzle. Thus, the initial 
and final shapes -will he similar, all dimensions being changed hy the 
ratio K of equation (l) . In most applications, the condition that 
the flow he -uniform and axiai. at the inlet and exit is desirable. In 
the region he-tween the -two ends, the cross-section shape -will generally 
•vary, althou^ this -variation has been found small in the flows computed 
up to the present time. It may he expected that the -variation of shape 
increases -with the rapidi-fcy of Mach number change between the inlet and 
exit. 


EXAMPLES 

Two axLsymmetrlc nozzles were designed hy the method of character- 
istics (see, for example, reference 3). The solution -was started hy a 
series expansion near the throat, -using the method of reference 4. 
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These nozzles -vrere designed xo expand from sonic velocity to a Mach 
number of 2.98 and 5 . 10 , respectively, and to have uniform axial inlet 
and exit flow. From these calculations, the resulting streamlines were 
found hy integrating along a characteristic to find the stream function. 
G3ae details of this computation are described in the appendix. In fig- 
ures 1 and 2 are given the resrolts of this calculation. 'With the use 
of these results, nozzles or inlets of any shape are readily formed. 

The details for the case of the Ifech mmiber 2.98 nozzle of figure 3 will 
be given. The exit section was selected as a square centered on the 
axis and having a side of length 14.5. The stream function was then 
found (at the exit) from equation (6) of the appendix for the periphery 
of the square. Throughout the len^h of the nozzle on the wall, ifr will 
be constant for constant values of the cylindrical angle. Thus, for 
the example considered, at the corner of the section, ijr = 13.10. Then, 
from figure 1 at x = 30, for instance, for this value of i{r, r is 
equal to 16.20. In this manner the shape of the nozzle may be found. 

The result is shown in figure 3. This illustrates how little the 
cross section varies from the initial shape for this particular basic 
flow. In figure 4, there is shown a similar nozzle designed for 
M =« 5.70. In this case it is seen that the nozzle shape shows approxi- 
mately the same variations as in llie case of the lower Mach number. 

Of course, if exit cross sections with only two flat walls are required, 
the other tvro could have circular-arc cross sections. Such a configura- 
tion wOTild probably be easier to fabricate than the square nozzles 
already discussed but would still have flat side walls suitable for 
schlleren windows. 

If it is desired to have perfect flow at the exit, such as is 
reqxiired in a wind t unn el nozzle, and still to have constant cross- 
section shape, the method of reference 5 might be xised to modify the 
cross sections already found. This possibility was pointed out by 
Ferri at the February 1952 meeting of the American Physical Society. 

On the other hand, if it is not essential to have perfect flow at the 
exit, sufficient accuracy might be obtained by the use of constant cross- 
section shape through the duct. 

Another iise of the method would be for the design of diffuser inlets. 
Such a configuration is shown in figure 5. In this case the inlet and 
exit sections are circles tangent to the axis of the basic flow. Since 
it was designed to diffuse to M 1, the inlet would have to be per- 
forated in 03:tler for the shock to be swallowed, but not as much as in 
the case of the usiial axisymmetric perforated inlet, because some 
external compression is obtained by cutting the diffuser back along 
the ^feich cone A (fig. 5). For the same reason, if the inlet were 
designed for only the Ksntrowitz contraction ratio (reference 6) so 
that no perforations were needed, the permissible contraetion ratio 
based on the free-stream tube swallowed would be increased. 
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To Illustrate this point, rou^ computations were made for two 
convergent-divergent diffusers, one symmetric and the other not. A 
conventional convergent- divergent diffuser designed for a Mach number 
of 2.98 has a permissible contraction ratio of 1.39 and a pressure 
recovery of 44.9 percent when the shock occurs at the throat Mach number 
of 2.63. The comparable diffuser designed by the present method would 
be one which appears as in figure 5, except that the area ratio from 
section B to the throat would be the Eantrowitz ratio corresponding to 
the average Mach number at section B (fig. 5) . To get some numbers, 
this average Mach number may be assumed to be 2.70. The permissible 
contraction ratio for this Mach number is 1.35 giving an over-all 
contraction ratio for the diffuser of 1.76 and a pressure recovery of 
54.8 percent when the shock occirrs at the throat Mach mmiber of 2.38. 

These two inlets have the same mass flow and entrance area, but the 
unsymmetrlc one has about 10 percent hi^er theoretical pressure recovery 
than the uaual convergent-divergent diffuser. 

The configuration shown in figure 5 could also represent a convergent- 
divergent exit jet. Such a jet might be useful in fairing the exit to 
the tail structure or aft portion of an airplane. 


CONCLUDING RS4AEKS 

A method has been presented for obtaining theoretical wall contours 
of supersonic nozzles and Inlets of arbitrary cross-section shape from 
knom axi symmetric or two-dimensional flows. If the basic flbw is 
axisymmetric and uniform at the throat and the exit, the desired nozzle 
will have similar cross sections at these stations. For the examples 
computed, the wall contoiir crass sections between the throat and exit 
did not deviate drastically from similarity with the initial shape. 


Le'sdLs Flight Eropulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, February 5, 1952 
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APEEKDIX 

DEOIEEMniATION OF STREAMLIHES AXIfflCMMETRIC FLOWS FROM 
KPTOWN VELOCITIES OR A CBARACTERISTIC HET 
The continuity equation in aod.al ly symmetric flow may te written as 


^ (p¥r cos ^ (p¥r sin 0) « 0 

where p is the cLensity, ¥ the speed, and the remaining quantities 
are defined by the sketch. 



Hence, there exists a stream function, ilr(x,r), such that 

^ ex p¥r sin 6 ^ ^ p¥r eos 0 

Sx pQ PQ 

where Pq is the stagnation density. Then, along the T)-characteristic 


Bt] 3x 3t] Sr St] 


or, from the sketch 

^ ^ cos (p + 0) + sin (H + 0) 

St] Sx or 


( 3 ) 


If equations (2) are substi't^ted into (3), there resvilts, after some 
manipulation. 


S^jf ^ p ¥r 
Sq Pq M 


(^) 


From the sketch. 


dq 


dr 

sin (p + 0) 
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Therefore, 




tB - 


gW 

PqM sin (p + e) 


r dr 


If the interval AB is small, say the distance between two points 
on the characteristic net, then it may he assumed, that the quantity 
oW 

1— takes its average. Experience shows that this will he 

PqM sin (p + 0) 

a particularly good, approximation along characteristics proceeding 
downstream away from the axis. Then 




B 


tA = 


7 \ I 

1 pW \“| 

\4pqM sin (p + e)J^ 

\4pqM sin (p + 0)/j 


^B 


- r 


L^) (5) 


At the exit, 0 « 0 and. M » constant. In this case, from the second 
of equations (2}, if the axis is taken as the i{r = 0 streamline. 


tB = 




( 6 ) 


Equations (5) and (6) are as accurate as the characteristic diagram. 
In the applications of this report, ¥ has heen made nondlmensional 
hy dividing hy the maximum speed. 
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